NPS69=82-005 


NAVAL POSTGRADUATE SCHOOL 


Monterey, California 


PERFORMANCE RATING OF ENHANCED 


MARINE CONDENSERS 


by 


R.H. Nunn 
P.J. Marto 


August, 1982 
Report for period ending 30 September 1982 


Approved for public release; distribution unlimited. 


Prepared for: 


David Taylor Naval Ship R&D Center 
Annapolis, MD 21402 


NAVAL POSTGRADUATE SCHOOL 
Monterey, California 


J. J. Ekelund, RADM, USN D. A. Schrady 
Superintendent Acting Provost 


This report documents progress for the period ending 
30 September, 1982, in the computer modelling portion of the 
project titled "Condenser Heat Transfer Augmentation." 

The work reported herein has been monitored and supported 
by the David Taylor Naval Ship Research and Development Center, 
Code 2721, Annapolis, MD 21402. 

Reproduction of all or part of this report is authorized. 


This report was prepared by: 


Wiililam M. tOLstes 
Dean of Research 


= Be A 


UNCLASSIFIED 
ES 
SECURITY CLASSIFICATION OF THIS PAGE (When Data Fntsrad) 


Teint mer rcen ‘ Lac READ INSTRUCTIONS 


; 1. REPORT NUMUER 2. GOVT ACCESSION NO 3. RECIPIENT'S CATALOG NUMULR 
NPS69-82-005 


4. TITLE (and Subittie) S. TYPE OF REPORT & PERIOD COVERED 
Progress report for 


COMPUTER RATING OF ENHANCED MARINE period ending 30 Sep 82 
CONDENSERS 6. PERFORMING ORG. REPORT NUMGER ; 
a 


iF. AUTHOR(s) 8. CONTRACT OR GRAWT NUMGER(a) 


R. H. Nunn | 
Pew. Matte | 
| 


9. PEHFORMING ORGANIZATION NAME ANO ADDRESS 10. PROGRAM CLEMENT, PROJECT, TASK 
AREA & WORK UNIT HUMBERS 
Naval Postgraduate School 


62543N; SF43400391 
Monterey, CA 93940 NO00167-82-WR-0114 


yil. CONTROLLING OFFICE NAME AWD ACDRESS 
DTNSRDC, Code 2721 
} Annapolis, MD 21402 


12, REPORT DATE 


August, 1982 


13. HUMBER OF PAGES 


92 


'S. SECURITY CLASS. (of thie report) ‘ 


WIONITORING AGENCY NAME & AODNESS(II different from Controlling Office) | 


UNCLASSIFIED | 


\Sa. DECLASSIFICATION/ DOWHGRACING 


SCHEDULE 


16. DISTRIBUTION STATEMENT (of thie iieport) 


Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (of tho obdotract ontered in Glock 20, If different from Report) 


1€. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continus cn reveree aide {{ neceseary and identify by block number) 


Heat transfer, Condenser, Enhancement, Computer modelling 


20. ABSTRACT (Continue on reverse side if neceesary oni icantify by Block manriser) 
Condenser tubes with enhanced surfaces show promise for 
use in advanced marine condenser designs. In order to ade- 
quately predict the performance of enhanced tubes when operating 
in a condenser-bundle environment, considerable improvement in 


| 
janalytical method is necessary beyond the standard approaches 
; 


+ recommended for use with conventional designs. 


| _ This report describes a computer code of moderate complexity 


om. FORM "as 
DO ssan yy 280g  EDIVION OF I WoY Ob IS GusoLE TE UNCLASSIFIED 
ee 


S/M 0102-°014* 4601 | 


SL akall SECURITY CLASSIFICATION GF THIS PAGE (Ghen Dale entered) 


UNCLASSIFIED 


a A ES 
Suc umly CL ass FICATION OF THIS PAGE(When Nore Entered. 


that is suitable for predicting the relative merits of the use 

of enhanced tubes in radial-flow condensers. The code accounts 
for vapor shear effects, condensate inundation, and noncondensi- 
ble gases, and predicts the local steam-side conditions at inter- 
mediate locations within the condenser bundle. 

Several computer experiments are described in which the 
impacts of enhanced-tube designs on overall condenser performance 
are investigated. The results serve to illustrate the importance 
of internal detail in such analyses and, in addition, they show 

‘that judicious use of enhanced tubes can lead to considerable 
improvements in condenser performance as well as the realization 
of heretofore unattainable design options. 


n 1473 UNCLASSIFIED 
4 7 Ss snes 
S/N 0102-914-6601 iv SECURITY CLASSIFICATION OF THIS PAGE(WHEn Date Entered) 


SUMMARY 


Condenser tubes with enhanced surfaces show promise for 
use in advanced marine condenser designs. In order to ade=- 
quately predict the performance of enhanced tubes when operating 
in a condenser-bundle environment, considerable improvement in 
analytical method is necessary beyond the standard approaches 
recommended for use with conventional designs. 

This report describes a computer code of moderate complexity 
that is suitable for predicting the relative merits of the use 
of enhanced tubes in radial-flow condensers. The code accounts 
for vapor shear effects, condensate inundation, and noncondensi- 
ble gases, and predicts the local steam-side conditions at inter- 
mediate locations within the condenses bundle. 

Several computer experiments are described in which the 
impacts of enhanced-tube designs on overall condenser performance 
are investigated. The results serve to illustrate the importance 
of internal detail in such analyses and, in addition, they show 
that judicious use of enhanced tubes can lead to considerable 
improvements in condenser performance as well as the realization 


of heretofore unattainable design options. 
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PERFORMANCE RATING OF ENHANCED 


MARINE CONDENSERS 


INTRODUCTION 


Since the early days of naval steam propulsion, the sur- 
face condenser has evolved into a very reliable component of 
the steam propulsion machinery plant. This reliability has 
been achieved by providing generous design margins to insure 
thermal performance at full power. The penalties for this 
overdesign, however, are additional weight and volume (which 
must be carried around for the life of the ship) and crowded 
machinery rooms with poor accessibility. Although surface 
vessels and submarines have different design constraints, both 
types of vessels can benefit from a more compact condenser de= 
sign through improved heat transfer. The actual dimensions 
of the condenser may have an impact on vessel size and dis- 
placement, the location of the propulsion machinery aboard the 
vessel, and even on vessel performance and cost. In addition, 
the future development of compact steam systems, both for main 
propulsion and secondary heat recovery purposes, will require 
the application of advanced technology to all system compo- 
nents, including the condenser. 

At the Naval Postgraduate School, the relative benefits 
of such advanced condenser designs have been analyzed using 


a modified, radial flow computer code (MORCON) to model the 


neat transfer and fluid flow phenomena which occur within 
this type of heat exchanger. Refinements to this code and 
improvements in the models used to predict two phase heat 
transfer and pressure drop are continuing. 

Based upon the work accomplished to date, it is clear 
that present day smooth=-tube steam condensers, when operating 
under typical conditions, are limited in thermal performance 
due to a large thermal resistance which occurs on the sea- 
water side (tube side) of the condenser. This resistance, on 
the average, is larger than the thermal resistances which oc- 
cur on the steam side and in the tube wall, and those due to 
seawater fouling and noncondensible gas buildup in the vicinity 
of the tube. Consequently, in present-day conventional con- 
denser designs, what happens on the steam side does not have 
to be modeled too accurately to arrive at reasonable designs.* 
When using enhanced tubes, however, the enhancement on the sea-= 
water sie can increase the inside heat transfer coefficient by 
100 to 200 percent over the smooth-tube case. The outside heat 
transfer coefficient, on the other hand, is increased by only 
LO to 50 percent. In this situation, the thermal resistances 
On the inside and outside of the tube can be approximately 
equal. As a consequence, it becomes essential to understand 
and accurately model the detailed mechanisms which occur when 


steam condenses on these enhanced surfaces. 


ate should be noted that, even with conventional designs, there 


can be locations within a condenser bundle where the steam side 
resistance is largest due to low vapor velocities and non- 
condensible gas pockets. 


Heat transfer on the steam side of Naval condensers 
often occurs in the presence of high vapor velocities and 
with severe condensate inundation from neighboring tubes in 
the bundle. It is also important to predict the resulting 
two phase pressure drop in order to trace the steam flow 
path through the condenser bundle and minimize losses in the 
saturation pressure and temperature across the bundle. Since 
comprehensive test data on this subject are virtually non- 
existent, and since advanced computer codes will depend 
strongly upon good steam-side models, additional experimental 
work is necessary. 

In spite of the tentative nature of existing theoretical 
steam-side heat transfer models, it has proved to be useful 
to conduct computer-based condenser performance calculations 
in order to provide comparative information regarding the rela- 
tive costs and benefits of enhanced-tube condensers. In addi- 
tion, the present theoretical basis is adequate to permit the 
identification of condenser tube-bundle regions that should be 
particularly sensitive to the introduction of internal and/or 
external enhancement. 

This report documents these computer-based calculations. 
Since steam-side conditions are especially important in such 
calculations, some preliminary discussion is devoted to this 
subject. This is followed by a brief general description of 


computer modelling of surface condensers. The computer code 


used at NPS is then described and several case=-studies are 
provided to illustrate the use of the code and to give some 
insights concerning the impact of enhanced-tube technology 


on condenser design. 
EFFECTS OF LOCAL STEAM-SIDE CONDITIONS 


Single-tube theories and experiments have shown, not 
surprisingly, that the thermal performance of a condensing 
tube may be dramatically affected by the thermo-fluid state 
of the surrounding two-phase multi-component flow. Thus the 
condensate loading and velocity of the vapor, its temperature, 
and its non-condensible gas content all have important influ- 
ences on the thermal resistance on the steam side of the tube. 
These influences lead to further complications in the analy- 
sis of condenser tubes performing in a bundle environment and 
the overall bundle performance is strongly coupled to the 
bundle inlet conditions, individual tube performance, and 
the general nature of the tube distribution within the bundle. 
Several of the main factors that are likely to be significant 
in bundle performance are discussed in the following para- 


graphs. 
LOCAL STEAM=SIDE PRESSURE VARIATIONS 


Frictional pressure losses on the steam side of a con- 


denser bundle are detrimental from a system point of view 


because these losses are directly chargeable against the 
thermal efficiency of the power system. In addition, the 
vapor temperature, which is a major factor in the heat trans- 
fer process that accompanies condensation, is strongly de- 
pendent upon the pressure in the vapor. This dependency is 
accentuated at low condensing pressures, varying from about 
7°F/in. Hg at 6 in. Hg condensing pressure to 12°F /in. Hg at 

3 in. Ho. (in other words, at about 3 psia a loss of 0.1 psia 
in pressure will lead to a 1.3°F loss in vapor temperature and 
at 1.5 psia the loss would be about 7a at ep Consequently, 
when condenser tubes are closely packed in a bundle the re=- 
sulting high vapor velocities and correspondingly high pres- 


sure losses can severely limit the bundle performance. 
VAPOR VELOCITY EFFECTS ON STEAM-SIDE THERMAL RESISTANCE 


High vapor velocities can have at least one effect that 
is beneficial. This is the shearing force that the vapor ex- 
erts on the condensate film which, at relatively high vapor 
velocities, can be significant. A number of investigators 
Cl, 2, 3, 4] have contributed to the understanding of the 
Serects of vapor shear. .In particular, Fujii, et al. [1] 
have recently proposed experimental correlations of the form 


CaSi= 2a) 


Re Keds) 


Nu, a fe) L 


where Nu, is the mean Nusselt number, Re. is the "two-phase" 


L 
Reynolds number (based on vapor velocity, tube outside di- 
ameter, and kinematic viscosity of the condensate), and Nu, 


is the standard Nusselt number for the ideal (zero shear) 


case. The empirical constants Cy and a are in the ranges 


gS Be Cy < 1.24 


kao Sa ee 

depending upon whether the tube thermal conditions are most 
appropriately described by constant temperature or constant 
heat flux conditions. The larger values are for the constant 
heat flux case. According to Fujii, et al., Eq. (1) is valid 
in the range 3.3 > Re, /Nu, > 0.28. For smaller values of this 
parameter, they recommend the use of a slightly reduced value 
of the standard Nusselt number that is given by 


Nu, = 0.96 Nu, 


In a typical naval condenser, vapor velocities well above 
100 ft/s are not uncommon and two-phase Reynolds numbers on 
the order of 19° are typical. For a nominal Nusselt number 
value of Nu, = 400, the resulting vapor shear effect is given 
by Nu_/Nu, ~ 1.5. Thus, the vapor shear effect is clearly 
important in the range of operation of naval condensers, and, 
when enhancements of various kinds are applied, this effect 
may become a dominant factor in overall condenser performance 


predictions. 


Correlations of the form of Eq. (1) appear to be ade- 
quate for the prediction of shell-side film coefficients in 
horizontal and downward vapor flow without condensate inun- 
dation. For upward vapor flow, however, the effects of vapor 
shear on film thickness appear to be much larger, and at very 
high velocities the shell-side film coefficient, normally low 
due to poor condensate removal, can be significantly increased 
- even to the extent that the influence of gravity is overcome 
by vapor shear effects. 

From the point of view of enhancement, simple reasoning 
would lead to the expectation that improvements in heat trans- 
fer due to vapor shear will be less significant at high ex- 
ternal enhancement conditions. This follows from the argument 
that enhanced tubes will carry thinner faster-moving conden- 
sate films (other conditions being equal) that are less 
sensitive to shearing effects acting at the vapor interface. 
Equation (1) reflects this effect in that (since a < .25 in 
all known cases) the ratio Nu /Nu, varies inversely with Nu,- 
The vapor shear effect might also be reduced by some of the 
previously-mentioned external enhancement schemes employing 
extended surfaces that would trap or otherwise "protect" the 
condensate. With high levels of internal enhancement, on the 
Other hand, failure to provide for adequate condensate removal 
will be somewhat compensated by the shearing of the relatively 


thick condensate film at high vapor velocities. 


Other vapor velocity effects, such as rippling, induce- 
ment of turbulence and stripping of the condensate film, are 
relatively poorly understood. A particularly important area 
in need of further research is the consideration of the en- 
trainment of condensate droplets in fast-moving vapor flows. 
As will be discussed in the subsequent section, condensate 
droplets as they traverse the condenser can have a profound 
effect upon the performance of downstream tubes. To the ex- 
tent that the paths of these droplets and the vapor flow are 
coupled through viscosity, an additional dependency upon vapor 


velocity must be considered. 
CONDENSATE INUNDATION 


The true behavior of the condensate generated within a 
tube bundle is far from the ideal model of vertically-dropping 
laminar sheets hypothesized by Nusselt. In fact, the conden- 
sate in most condenser environments will be in many forms 
(droplets, rivulets, sheets, sprays) and will flow in many 
directions. It is highly unlikely that any quantity of con- 
densate will traverse a condenser without losing its identity 
by impinging upon other tubes in the bundle. The effect of 
this impingement is referred to here as condensate "“inunda- 
tion" or "rain" and it is probable that a true theoretical 
model of the phenomenon would of necessity be stochastic in 


nature. In the absence of such analytical sophistication, a 


number of empirical models have been developed that are useful 
when applied within their range of credibility. 

Strict application of the classical Nusselt theory, to- 
gether with its restrictive assumption of a laminar film fall- 
ing under the influence of gravity alone, leads to the conclu- 
sion that the average condensate film heat transfer coefficient 


for a vertical row of n tubes, hoe is given by [5] 


i. 70-25 (2) 


where he is the single-tube film coefficient. Numerous inves- 
tigators have compared this simple expression with experimental 
results [6, 7, 8, 9, 10, 11] and it has become clear that in 
reality the results obtained from Eq. (2) reflect an over- 
prediction of the detrimental effects of inundation. In fact, 


a modified form of Eq. (2) may be written as 
h 
° 


where empirical evaluation of s has led to values in the range 
“yes <'525 [7]. Within this range, for example, a column of 
30 tubes may be judged to have an average film coefficient that 
is 79% to 43% of the single-tube value. It should be noted 
further that the actual film coefficient for the n-th tube in 


a column, h,, may be found using Eq. (3) with the result that 


ae 3 n{i-s) 


Ag 


(l-s) 


(n-1) (4) 


The range of s leads to a range in ho of from 75% to 32% of 
ho: clearly the range of uncertainty ins is unacceptable. 

Among the chief discrepancies in the assumed constancy 
in s in Eqs. (3) and (4) are: 

l. The effect of inundation is definitely dependent 

upon the vapor velocity and direction. 
2. Bundle tube layout (tube spacing and pattern) will 
influence condensate inundation. 

Thus, even if the functional relationship given in Eq. (3) is 
authentic, the exponent s must show the dependency 

s = fct (steam velocity vector, bundle geometry). 
The first of these is indicated nicely by Brickell [11] in 
which a general decrease in s with steam velocity is shown. 
These results are reproduced here in Fig. 1. It is also ap= 
parent that at low steam velocities the exponent s is larger 
(more inundation effect) for both upward and downward steam 
flow than it is for horizontal steamflow. This directional 
effect is less noticeable at high steam velocities. The in- 
fluence of tube separation upon inundation effects is not 
clearly discernable from available data, largely because of 
the coupling of vapor shear effects with those due to inunda- 
tion. It is possible to speculate with some confidence, how- 


ever, that closely-spaced tubes are more influenced by 


inundation than are the tubes in relatively disperse bundles. 
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One important influence of tube spacing in low velocity 
steam flow is that of side-drainage whereby condensate gen- 
erated from above may, due to surface tension effects, pro- 
ceeds laterally to adjacent tubes rather than directly down- 
ward. Eissenberg [7] has contributed a means of quantifying 
this effect with the result that Eq. (3) may be modified as 


follows: 


= tle Panne 


+ FAL006 + 2/2 (2n) 7] (5) 
where Fy depends upon tube spacing and approaches unity for 
densely—packed staggered tube banks and zero for disperse tube 
layouts. Equation (5) represents the state-of-the-art in the 
prediction of condensate inundation effects and its effective- 
ness in doing so is, of course, highly dependent upon the 
selection of correct values for s and Fa: 

The extent to which enhanced tubes are susceptible or 
resistant to inundation effects is a matter that is far from 
settled. (Indeed, the issue is unresolved as far as smooth 
tubes are concerned, when bundle effects are taken into ac- 
count.) There are some experimental results [12] that indi- 
cate that smooth tubes operating in the dropwise mode of 
condensation are not greatly affected by inundation, (This 
observation was made at very low vapor velocities.) Prelimi- 
nary indications are likewise encouraging for film condensation 
on enhanced tubes where single-tube experiments [13] have indi- 


cated values of s [Eq. (3)] in the range of 0:04 to 0:10; 
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NON=-CONDENSIBLE GASES 


The effects of non-condensible gases can be classified 
into one of two categories: 

1. The introduction of an additional local thermal 

resistance due to the propagation of gases towards 
a condensing tube surface under the influence of a 
gradient in partial pressure, and 

2. The cumulative effect of gas blanketing where uneven 

rates of condensation in a condenser bundle eventu- 
ally lead to regions where tubes are inoperative in 
a condensing role. 

In naval condensers [14] as well as in shore-based power 
condensers [8], non-condensible gases contained in the inlet 
flow are usually on the order of 100-200 weight-parts-per= 
Million. This ratio increases in the condenser, as steam is 
extracted, but even at a worst-case condition (20,000 wppm, 
say) the local overall heat transfer coefficient will be re- 
duced by only about 15% and this in a very small region of 
the condenser. Thus, although a number of theories are avail- 
able to account for the local effects of non-condensible gases 
[15] in naval condensers, there appears to be little need for 
extreme precision in the calculation. 

The effect of non=-condensible gases on the external re- 
Sistance of enhanced tubes is a matter of considerable uncer- 


tainty. Cunningham [16] has recently reported on single-tube 


1S 


ee in which fine-finned tubes were found to be less 
affected by non-condensible gases than were smooth tubes 

under corresponding conditions. In any event, it is ex- 

pected that because of the low concentrations of non-condensible 
gases under naval condenser conditions their effect will be 
small unless their presence in some way affects the enhancement 
process itself. 

The more-important aspect of non-condensible gas effects, 
crucial in all condenser applications, is the design of tube 
bundles such that non-condensible gases do not accumulate. 

In large tube bundles it is often necessary to incorporate 
special inlet baffling and steam lanes that allow the steam 
to penetrate deep into all parts of the tube nest in order to 
avoid gas blanketing. The compensation of condenser bundle 
designs for this purpose has often been experimental in na- 
ture, involving large-scale and extremely expensive tests, 
and is probably a major reason why condenser cross-sections 
do not seem to change drastically from generation to genera- 
ELON. 

In recent years, however, the ascendancy of computational 
capabilities has led to the feasibility of the comprehensive 
computer codes necessary for the prediction and alleviation 
of gas blanketing. One such computer code has been developed 
in the UK (such codes are not, in our experience, available 
to the public) and has been used to great advantage in the de- 


sign and improvement of shore-based power condensers [8]. 
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SUMMARY 


The performance of a condenser tube within a tube bundle 
is strongly influenced by the history of the vapor/gas mixture 
prior to arrival at the tube in question. Conversely, the 
performance of any given tube can have a profound effect upon 
the thermal behavior of tubes located downstream. The com- 
plexity of these interactions may lead to a degree of depres- 
sion in the mind of the researcher who would predict overall 
bundle performance. On the other hand, these inadequacies in 
understanding may be viewed as challenges and, even in our 
present state of ignorance, we can answer questions ina 
relative way. For instance, qualitative trends due to aug- 
Mentation may be identified and relative improvements due to 
augmentation may be quantified within a precision of perhaps 
10-208. | 

Just as a condenser is larger and more complex than a 
single tube, it is to be expected that the necessary compu- 
tational and experimental scales will be larger for condensers 
than they are for single tubes. To justify such investments, 
and to guide the efficient allocation of resources, it is 
necessary to use what tools we have at hand to investigate 
potential payoffs and pacing problems. The following sections 
describe some preliminary results obtained from computer 


models of condenser bundles. 
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COMPUTER MODELING 


Analytical models for the estimation of condenser per- 

formance may be classed as follows: 

1. Standard methods of the "black-box" type that are 
used for preliminary - and sometimes, it is sus- 
pected for final - design purposes. These methods 
include those recommended by the Heat Exchange In- 
stitute (HEI [17]) and British Electrical and Allied 
Manufacturers' Association (BEAMA [18]). 

2. Intermediate computer codes that, to varying extents, 
include the details of tube-to-tube flow within the 
condenser. For purposes of classification, these 
codes are restricted to the analysis of condenser 
designs in which the vapor/gas mixture may be pre- 
sumed to flow in a fixed direction (such as radial, 
horizontal, upwards, or downwards). The code MORCON 
is of this class and is described in some detail be- 
low. 

3. Multi-dimensional codes utilizing finite element or 
finite difference methods for the solution of the 
governing conservation equations. Because codes in 
this category are proprietary (in all cases, to our 
enowledeess little is known concerning the extent to 
which they exist and work. The CERL code mentioned 


above [8] appears to be quite successful in application 
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and there are some indications [19] that U.S. manu- 

facturers possess a similar capability. 

The black-box methods are extremely limited to their 
application to all but the most standard condenser designs. 
In essence, only the overall heat transfer coefficient, U, is 
predicted by these methods and even this is subject to con- 
siderable uncertainty (HEI predictions of U have been ob- 
served to differ from measured values by from 100% on the 
high side to 50% on the low side [20]. These methods of pre- 
diction are based upon the assumption that the overall thermal 
resistance is dominated by the internal (tubeside) resistance 
- they are limited to cases of moderate cooling water velo- 
city - and the allowances for fouling are subject to consid- 
erable criticism on physical bases. No means are provided 
for accounting for the distribution of thermal resistances. 
across individual tube sections or between tubes in a bundle. 
It is clear that major advances in condenser designs cannot 
be based on these black-box methods. 

At the other extreme of the spectrum of sophistication 
lie the multi-dimensional computer-based methods. The basic 
computational difficulties in the development of these com- 
prehensive computer codes are largely ones of scale. It is 
quite difficult, yet most important to include two-dimensional 
capabilities that allow calculations in the plane of flow 


normal to the axis of the condenser. The third dimension, 
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steals flow in the direction of the tube axis, may be ob- 
tained with reasonable accuracy by combining two-dimensional 
slices. Even so, the complexity of the two-dimensional cal- 
culation is sufficient unto itself in that it is often neces- 
sary to compute the two-phase multi-component flow in the 
passages between thousands of tubes. As a result, it is 
common to combine groups of tubes into "mini-bundles". Within 
these mini-bundles the equations of mass, momentum, and energy 
are satisfied by means of finite-difference expressions for 
the convective, conductive, and source effects. The source 
terms in these expressions require much of the approximation 
necessary in less-complicated codes in terms of interaction 
effects such as vapor shear and condensate inundation. Thus, 
the extension to multi-dimensional computations does not ob- 
viate the need for additional insights into the various bundle 
interactions mentioned in the previous section. The chief 
benefit of such extensions appears to be the ability to pre- 
dict the performance of poorly designed condenser bundles 
especially with respect to air-blanketing. 

Multi-dimensional computer codes can be extremely valuable 
for condenser performance predictions and may even be essential 
to the development of truly optimized condensers. Large re- 
source investments are required for the development of such 
codes and questions inevitably arise concerning their cost- 
effectiveness. This issue is not developed further here except 


to say that, because of it, our emphasis has been at the 
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intermediate level of analytical sophistication where codes 


such as MORCON are found. 
THE MORCON COMPUTER CODE 


MORCON (Modified Oak Ridge Condenser) is a modified ver- 
sion of the ORCON1 code developed at the Oak Ridge National 
Laboratories [21]. Some of the features and limitations of 
MORCON are summarized below: 

a. MORCON allows the calculation of radial flow in cir- 
cular bundles. (Other unmodified versions of ORCON1 
have been developed for downward and horizontal-flows.) 

b. Correction of the condensate film coefficient for the 
effects of rain, including side drainage, is accom- 
plished by the means previously described. As has 
been pointed out, the correction only accounts for 
the distribution of condensate rain upon lower tubes 
in a gravity-dominated flow and does not account for 
the interaction of vapor shear effects with those 
due to inundation. 

c. Non-condensible gases (up to four species) are taken 
into account by the use of the Colburn analogy [22] 
together with an empirical correlation credited to 
Eissenberg in [21]. 

d. MORCON allows for the effects of heat transfer en- 
hancement of the inner and/or outer film coefficients. 
The influence of enhancement on tubeside and shell-side 


pressure drops is also included. 
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MORCON, as in all other known cases, includes foul- 
ing effects only as a constant resistance term. 

All tube-side and shell-side fluid properties are 
continuously calculated as they change with pres- 
sure and temperature throughout the condenser. The 
latter quantities vary due to inlet, outlet, and 
tube-to-tube flow cross-section variations, as well 
as viscous losses. 

Bundles can be segmented into as many as six compu- 
tational sectors. Steam flow entering each sector 
is adjusted until the shell-side pressures at the 
exits of each sector (which is also the air cooler 
inlet) are matched within a specified limit. The 
sector arrangement for the circular bundle with 
radial flow is illustrated in Fig. 2. 

The overall heat transfer coefficient for each row 
in each sector is calculated on the basis of an 
"average" tube. This is a basic limitation of the 
code that is somewhat relieved by the ability to 
calculate performance by sectors. 

Vapor shear effects are included in the manner pre- 
viously described. To the extent that the Fujii 
correlation does not include the interaction of 
vapor velocity with condensate inundation, the re- 
sults are suspect for other than gravity-dominated 


(low vapor velocity) flow. 


20. 


air removal point 


Sec LOF 


Figure 2. Sketches showing 
sector arrangements and 
quantities necessary for 
inundation calculations 
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DESCRIPTIVE SUMMARY 
In this section the organization and flow of MORCON are 

described. A listing of the code is provided in Appendix A. 
Because MORCON is in large part identical to the ORCONI code, 
the treatment here will be brief with details confined to the 
few modifications made to ORCON 1. For a more complete des- 
cription of the basic code, the reader is referred to [21]. The 
subroutines contained in the code are listed below: 

MAIN (not labeled) 

SECALC 

COOLEX 

HETTRN 

ADJUST 

PRSRDP 

DFSVTY 

INPUT 

OUTP1 

ZERO 
Subroutine DFSVTY computes the mutual diffusivity of the vapor/ 
gas mixture. Various function subroutines are included to cal- 
culate specific heats, thermal conductivities, viscosities, and 
densities of the gas, vapor, and condensate as required. Func- 
tion subroutines are also provided to calculate vapor saturation 
pressures and temperatures and the latent heat of vaporization. 


The purposes of subroutines INPUT, OUTP1, and ZERO are indicated 
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by their names and these will receive no further discussion here. 


The subroutine calling protocol is indicated in Table I, below. 
fabie: i... Subroutine Calling Protocol 


Calling Called Subroutine 
Sub= 
Routine j| SECALC | COOLEX | HETTRN | ADJUST | DFSVIY | PRSDRP INPUT |OUTPL 


var is i AS SS 


Subroutine MAIN 

Subroutine MAIN is the main control program. Input values 
are taken into this subroutine and the tube bundle cross-sectional 
geometry is determined. This geometry is defined by the specified 
number of tubes, their outside diameter, and the spacing between 
them. The tubes are arranged radially such that the outer row of 
the bundle is complete. The tube array is staggered with the 
radial spacing between rows equal to (/3/2) times the specified 
circumferential spacing between tubes. The void diameter is de- 
termined by the space necessary to accommodate the number of tubes 
that are insufficient to complete an outer row. With the tube 
sheet entirely filled (zero void diameter), the specification of 
additional tubes will result in the removal of a number of tubes 
from the void sufficient to complete an additional outer row. 


(Because of the difference in performance of inner and outer-row 
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tubes, this method of allocating tube cites leads to highly 
non-linear and sometimes confusing results when tube count is 
used as an independent design variable.) 

Bundle symmetry is assumed about the vertical centerline 
of the circular bundle and the remaining semi-circular tube 
array is divided into six sectors. Each partial row in each 
sector is defined by the central tube in that string and the 
program calculates the number of tubes in a vertical row above 
this "mean" tube. This is needed to provide for inundation 
calculations. 

Following these preliminary arrangements, subroutine 
SECALC is called. The return information is the condition of 
the vapor/gas mixture as it leaves the main condenser bundle 
and enters the air cooler. The performance of the cooler is 
then computed by means of a call to subroutine COOLEX. 

The final overall vapor/gas conditions, as returned from 
COOLEX, are transmitted to subroutine ADJUST. If the exit 
Fraction of the current design is equal to the input specifi- 
cation (within five percent error) the return from ADJUST is 
followed by final output actions. If the specified exit frac- 
tion is not obtained, the code will adjust either the inlet 
steam flow or the tube (bundle) length, as selected by the 
user. The updated design is then returned to SECALC for evalu- 
ation and the above procedure is repeated until convergence is 
obtained or the iteration limit (set to 10) is exceeded. 


Figure 3, below, summarizes the function of subroutine MAIN. 
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MAIN 


CALL 
INPUT 


DEFINE 


GEOMETRY 


CALL Ps 
CECAILG condenser outlet conditions 
GAEL ~~ 
COOLEX cooler outlet conditions 
CALL change steam flow or 
ADJUST tube length as required 
ES CALL 
OUTP1 


Figure 3. Simplified flow diagram of main calling subroutine. 
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Subroutine SECALC 

Subroutine SECALC, with calls to subroutine HETTRN, de- 
termines the sector-by-sector and row-by-row performance of 
the condenser. Calculation is begun in the first (uppermost) 
sector at the first (outermost) row. Using the central tube 
of this and subsequent row segments, subroutine HETTRN is 
called to determine the thermal performance of the current 
"man" tube. With information returned from HETTRN, subroutine 
PRSDRP is called to determine the pressure of the vapor/gas 
mixture entering the next row and the steam quantity and con- 
ditions for this next row are returned to HETTRN. 

The outer loop in SECALC covers the six sectors while the 
inner loop covers the rows (the number of rows is determined 
by the input geometry). Following the completion of both 
loops, the steam=-side pressures at the outlet of each sector 
are each compared with the mean of the outlet pressures. The 
total inlet steam flow is then adjusted so that sectors with 
excessive steam-side pressure drop receive less flow and con- 
versely for those sectors with relatively higher outlet pres- 
sures. The calculation is then repeated in its entirety (all 
sectors, all rows) until the bundle outlet pressure is uniform 
to within one percent. 

Figure 4, below, illustrates the function of subroutine 


SECALC. 


SECALC ) (called by MAIN) 


UPDATE 
MIXTURE 
OND Chi 


CALL “mean” tube thermal performance 
HETTRN 

CALL shell-side pressure drop 

PRSDRP 


ADJUST : 
YES 
FLOW TO sete {RETURN 
Uk 


Figure 4. Simplified flow diagram for subroutine SECALC. 
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Subroutine COOLEX 

Subroutine COOLEX is a simplified version of subroutine 
SECALC in which calls to HETTRN are used to obtain row-by-row 
performance of the cooler tubes. No steam-side pressure 
balance (such as that performed in SECALC) is attempted in 
COOLEX. 

The user specifies the percentage of the total number of 
condensing tubes to be used in the cooler section. These 
tubes are allocated to a vertical staggered array with cross- 
flow dimension initially equal to the bundle void diameter. 
If this tube cite allocation results in crowding of the cooler 
tubes (as indicated by vapor/gas mixture velocities in excess 
of 150 ft/s) the cooler is widened (and proportionately 
shortened) as necessary. 

Subroutine HETTRN 

This program contains the basic heat transfer algorithms 
and, as noted above, is called by SECALC for each mean tube 
in each row of each sector. Following a summary of the pro- 
cedures followed in HETTRN, a few of the more-important heat 
transfer algorithms are discussed in some detail below. 

Because the thermodynamic and physical properties of the 
cooling water are allowed to depend upon temperature in HETTRN, 
a basic iteration is required to compute the overall heat 
transfer coefficient. This iteration begins with an estimate 


from which, with local steam conditions and cooling water inlet 
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temperature and flow rate, the cooling water outlet tempera- 
ture is calculated. Cooling water properties are then calcu- 
lated based upon the mean of the inlet and outlet temperatures 
of the cooling water. 

Following the determination of cooling water properties, 
HETTRN then computes, in order, the internal thermal resistance 
(based upon theDittus-Boelter correlation with internal en- 
hancement taken into account) and the external condensate 
film thermal resistance based upon the Nusselt correlation 
(see below). The external resistance is then corrected for 
vapor shear and condensate rain effects. These resistances, 
together with tube conduction resistance and a constant 
fouling resistance (input by the user), are then summed to 
develop an interim value of total resistance to heat transfer, 
exclusive of the resistance of the non-condensible gas film. 
The thermal resistance of the non-condensible gas film is 
calculated and various key variables for the current iteration 
are determined. Convergence of the basic loop is based upon 
agreement between current and previous values of film tem- 
perature drops. Agreement to within 0.01 degrees Fahrenheit 
is required in the case of the liquid film temperature drop 
and if non-condensible gases are present the iterations con- 
tinue until the calculation of the gas film temperature drop 
is steady to within 0.001 degrees Fahrenheit. Following con- 
vergence, values are transferred to SECALC and common storage 


as necessary for future use. 
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Condensate Film Thermal Resistance 
In MORCON the standard Nusselt relationship for the 
condensate film resistance is corrected for vapor shear ef- 
fects according to the correlations of Fujii et al. [1].° Gn 
the assumption that the thermal resistance on the steam side 
is relatively small, the constant heat flux case is taken to 
be most appropriate. (Within experimental uncertainties the 
two cases = constant heat flux or uniform wall temperature - 


are not appreciably different.) Eq. (1) becomes 


AS Oe Nui? Re, °* (6) 


where Nu, is the Nusselt relationship given by 
v 
cpm aisle : 
= 0.2 (ay 
Ku, AT 


Following correction of the mean Nusselt number for 
vapor shear effects, this value is further corrected for 
condensate inundation. Equation (5) is applied with a value 
of s = 0.2 so that the average film coefficient for a ver- 


tical row of m tubes is given by 


h (02 
= 0.6F 3 + (l - 0.5647F,)m (8) 


This expression is evaluated for m = n and m = n-l, where 


nis the number of condensing tubes vertically above the 
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tube in question, and the heat transfer coefficient for that 


tube is given by 


Dei Boe (n-d)ehe (9) 


n n 1B 


The selection of the exponent s = 0.2 in these calcula- 
tions is based upon the comparison with experiment shown in 
Fig. 5 and is also in consonance with the observations of 


eeuhers [7, 11, 23] (see also Fig. 1). 


Figure 5. Tubes with 
inundation and side 
drainage; comparison of 
theory with experiment. 
Data of Ferguson and 
Oakden ("Heat Transfer 
Coefficients for Water 
and Steam in a Surface 
Condenser," Trans. Chem. 
Eng. Congress (World 
Power Cont.), v.37; 640 
(1936)) as reported in 


Experiment 
O Tubes in bundle Ea 


O Column of tubes 
— Theory, Eq. (8) 


Non-Condensible Gas Film Coefficient 
The process of heat transfer through the non-condensible 
portion of the gas/vapor mixture is modelled as one due to 


mass transfer of the vapor through the gas. Sensible heat 
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transfer across the gas film (usually about 0.1% or less) is 


neglected. A mass transfer coefficient, K is defined such 


(ey 


that BasAT h 
K = ee ee ( 
ri 


— (10) 
G h g Pst Neg 5 | 


where AP. is the drop in vapor partial pressure in the 


= 
film, AT is the corresponding temperature change, Reg is the 
latent heat of vaporization, and he is the heat transfer co- 
efficient for the process across the non-condensible film. 

The drop in vapor partial pressure across the gas film 
is taken equal to the rise in gas partial pressure across 


the. fiim,so that 


Ae g joe — Pas 


(12) 
where Pc and P ob are the gas partial pressures at the con- 
densate film interface and the external mixture conditions, 
respectively. Combining Eqs. (10) and (11), the gas film 


temperature drop may be written: 


Ken 
Rey ie a=) - P oT 
AT a, gc gb (3 | (12) 


The partial pressure P ob is known at each point in the calcu- 
lation and the rate (dT/dP) is calculated from the saturated 
steam line. The condensate interface gas partial pressure 
is taken from the empirical relationship suggested by Eis- 


senberg and Noritake [24]: 


Bal 
Pp = 
gc ve 
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where 


ee eS G exp (.54 - .544 ln Re) 


/ 
(so) 3 


with G the local mass velocity, Sc the local Schmidt number, 


and Re the local vapor Reynolds number. With this correla- 


tion, Eq. (13) becomes 
Cn 
pape INES beh a 
AT hg Pg» (3 (13) 


The gas film temperature drop may be further defined in terms 


of heat transfer coefficients by the following heat balance: 


at = U AT 


g (14) 


where Uy is the overall heat transfer coefficient belonging 


to the log-mean temperature difference AT a aie: AsS) pao COs 


the heat transfer coefficient without non-condensible gases 


(known, at this point in the calculation) then: 


al iL GSES) 
Th Sey 


Combination of Eqs. (13) = (15) with the elimination of AT and 
Ag leads to the following quadratic relationship for the over- 


all heat transfer coefficient: 


2 hts. Cx P Ais Gold 
ts = (Pals sof ae (B)[}o, WE» Oe 
m 


2S) 


All the quantities necessary to determine the coefficients 
in this relationship are available (for the current iteration) 
soOwthat’ a solution for Us is forthcoming. Solutions for Ae and 
AT follow according to @gs..0(15) “and (14). 

Subroutine PRSRDP 

As the name implies, the purpose of this subroutine is 
to compute the pressure drop on the steam side between adjacent 
rows of tubes. In MORCON (as in ORCON1) the friction factor is 
given by [24]: 


£ = 0.102) toe oe (17) 
Re 


where £ is related to the pressure drop by 


4 as 
Gowte tale 
(ayig 
AP = 2£G? 
fe) 


The application of the correlation of Eq. (17) is limitediee 
moderate mass velocities (G < 1500 lb/hr - ft?) and, as such, 
may lead to significant errors when used to predict steam=-side 


pressure drops in some advanced marine condenser designs. 


SAMPLE RUN 


In this section a sample run is described in order to il- 
lustrate the necessary user procedure and to provide a base=- 
line case for subsequent analysis of the effects of enhancement. 


The tubesheet layout for the baseline design is shown in Fig. 6. 
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Figure 6. Assembly and tubesheet layout for baseline 


condenser. 
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and general specifications for this design are listed in 
Tablerli< 


Table II. Design and Performance Specifications for 
the Baseline Condenser 


Inlet steam: 189,003 lb /hr saturated 5 in.Hg (133.7°F) 
so.7 lb/hr noncondensibles (air) 

Tubes: 1646 tubes in each of two passes, 6% in cooler 
70-30 Cu Ni, 5/8 .OD¥ 16 -BWG (0 .06Sewadas 


Pitch/OD 1.35 4 2% 


Fouling Resistance 3. 32800 7 HELEC F/Btu 


Cooling water: Inlet temperature 66.1°F 


Velocity 8 ft/s (7,900 gpm) 
Salt concentration 3.5% by weight 


User Supplied Data 


Subroutine INPUT reads, writes, and processes the data 
rovided by the user. These data are supplied via seven 
Cards, as illustrated An’ Fig. 7. “The first five of these 
cards are identical to those described for ORCON1 [21] but 
the input nomenclature is repeated here so that the reader 


may proceed (it is hoped) without reference to that document. 
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Following is a list of definitions of the input para- 


meters. A listing of the input cards for the baseline case 


is given in Appendix B (compare Table II). 


Gard 1, (FORMAT 2. of1 ,«l9A4)0A3) 


IR 


INDENT 


Indicator for cards to be read (column 1 
only). 2f IR = 0, 4, or >5, then ail) seven 
cards will be read in proper order. 

All the remaining columns of the first card 
are available for identification or notes, 


79 columns 


Card 2, (F@RMAT: 7F10.0, 215) 


ANTP 


PRCCLR 


SDD 


SDDMIN 


SKW 


NRWSP 


Total number of tubes in the bundle (includes 
cooler). 

Percent of tubes in the cooler. 

Ratio of tube circumferential center=-to-center 
spacing to tube diameter 

Ratio of tube spacing to tube diameter for 
cooler section. 

Outside diameter of all tubes, in. 

Wall thickness of tubes, in. 

Thermal conductivity of tube material, 
Btu/heste-- one. 

Number of rows of condenser coolant inlet tem- 
peratures to be read (use only to read punched 
card inlet coolant temperatures, set = 0 other- 


wise). 
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NRWSCP Same as NRWSP except for cooler. 
Gard 3, (FORMAT: 8F10.3) 

WBI Total coolant flow to condenser (cooler in- 
cluded), lb/hr. 

VELBIP Coolant velocity, ft/sec. 
(Note: WBI and VELBIP are exclusive variables; 
one of them should be set to zero. If not, 
then VELBIP will override WBI). 

CBI Coolant salt concentration, fraction by weight. 

STBIP Coolant inlet temperature, so 
(Note: STBIP is used only if inlet coolant tem- 
peratures are not read in or transferred.) 

FOUL Tube fouling resistance hr-ft?-°F/Btu 

ENHI Tube internal heat transfer enhancement factor. 

Card4.,.9 (FORMAT s); SAL 7 PS 0:2),, 94:5) 

GAS Type of noncondensible gas, must be one of AIR, 
CO2),) fOr IMEX 

SECFLG Bundle model indicator, gives number of 30° sec- 
tors in one-half of a symmetrical bundle, or 
total sectors in a semicircular bundle. 

HFCDFL Symmetry indicator. 
If = 1, model will have symmetry about a ver- 
tical centerline 
If = 0, model will be a semicircular type. 

RADFLG Factor to increase bundle radius beyond that 
normally calculated by the code (will enlarge 


center void). 


Bye 


BAFFLE 


FDAVE 


EXITFR 


Q@UTPUT 


INSTM 


ITRAN 


IPLYAT 


Flag to indicate presence of baffles at 2 
and 4.0'clock. 
Tube spacing parameter. Used in subroutine 
HETTRN to modify effect of condensate flood- 
ing. 
Target value for vented steam as a percent of 
input., > Lf EXITPR= "0 .0;. then, program wii 
make a single pass and return the results 
without any adjustments. 
Control flag to restrict amount of detail 
printed out. 
GUTPUT = 1.0 gives full detail of all sectors. 
@UTPUT = Q.0 prints only summary. 
Flag to direct program flow when converging on 
exit fraction (EXITFR) target. 


INSTM 


1, inlet steam adjusted. 

INSTM = 0, tube length adjusted. 
Flag to direct code to use outlet coolant tem- 
peratures from immediate preceding run (used 
only when multiple stages are run together) 
for inlet temperatures of current stage. 

ITRAN = 1, use transferred temperatures 


ITRAN 0, do not use transferred temperatures. 


Flag to give option of having detail output in 
floating point form; one gives floating point, 


zero gives fixed point. 


40 


IPNCH Outlet coolant temperatures will be punched 
Porvainput, wt, LPNCH? =).0. 


Card 5, (FORMAT: 8F10.3) 


WSI Steam flow to condenser, lb/hr. 

WNCI Noncondensible gas flow to condenser, lb/hr. 
STSAT1 Inlet steam temperature, OF 

ALST Tube length, ft. 

ENH@ Tube external thermal enhancement factor for 


heat transfer. 
ENHF Friction factor enhancement factor for use in 
calculation of pressure drop in the steam. 
Card 6), “(FORMAT s). 4F10.3) 
ENHFI Friction factor enhancement factor for cooling 
water pressure drop 
TUBESW Specific weight of tube material, 1l Vee 
FC Coefficient for enhanced-tube friction factor 
formula of the form f = FCs (Re, /10>) **FEXP 
FEXP See FC 
Card 7, (FORMAT: 215) 
EVAP, Switches to control vapor shear and inundation 
IRAIN 
effects. If set less than one these effects 


are neglected. 


Output 


Subroutine OUTP1 handles the principal output from the 
run. First, an image of the input deck is reproduced followed 
by an expanded version of this information. There next follows 


a series of statements defining the progress of the iterations 
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to aentiats the specified exit fraction. In the baseline case 
(Appendix C) the tube length is successively adjusted until 
convergence is achieved in the ninth iteration. There fol- 
lows a table of data describing the condenser, cooler, and 
overall performance. Next, average inlet and outlet coolant 
temperatures are calculated for the condenser and the cooler. 
From these averages overall log mean temperature differences, 
DTCND2, DTCOL2, and DLTOT2, are calculated for the condenser, 
cooler and overall, respectively, using inlet steam tempera-= 
tures. For the appropriate values of heat removed, the overall 
heat transfer coefficients, UPCOND, UPCOOL, UPAVG, are found 
for the condenser, cooler, and entire bundle. In addition, 
during the run, a row=-by-row average of the heat transfer co- 
efficients (UBARW, UBARWC, UAVGW, respectively) is made, and 
the corresponding log mean temperature differences (ADTCND, 
ADTCLR, ADTOA, respectively) are back calculated. 

The option exists to print only the summary of results 
(OUTPUT = 0), or to have the summary plus two pages of detailed 
results for each sector (OUTPUT = 1). In the event of an error 
condition existing when the output subroutine is entered, the 
above calculations will be omitted and only the detailed output 
will be printed. 

Appendix C contains the summary output as well as the de- 
tailed output. In addition, .a listing of definitions /(moe 
normally included) is inserted to facilitate the analysis of 


the detailed output. At the bottom of the summary output 
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(preceding the detailed output) are seven additional quanti- 


ties. These are defined as follows: 


BUNL2P Tube length for two-pass bundle, ft 

BUNVOL Bundle volume, including tubes and voids, Be 
BUNWT Bundle weight, lb 

RETUB Cooling water Reynolds number 

FRIFAC Cooling water friction factor 

HDLOSS Cooling water head loss, ft H50 

PMPPWR Cooling water pumping power, hp 


RESULTS OF COMPUTER EXPERIMENTS 


The following calculations examine applications of the 
MORCON code to departures from the baseline design. Several 
computer experiments are described in order to demonstrate 
some of the effects of internal and external enhancement of 
tubes contained within bundles. In addition, these examples 
serve to illustrate some of the influences that tube=-to-tube 


interactions can have upon the performance of tubes in bundles. 


Effect of Internal and External Enhancement on Bundle 

Length 

Figures 8a and 8b illustrate the predicted effect of 
internal and external enhancement in terms of the length of 
tubes necessary to Maintain the baseline heat duty. Constant 
heat duty is obtained by iteratively varying the tube length 


EOMobtain a fixed exit traction of 0.05% of the baseline. inlet 
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steam flow. Figure 8 includes manufacturer's data for en- 
hancement factors of some commercially available augmented 
tubes and it is seen that combined internal and external en- 
hancement can be expected to lead to bundle length (and volume 
and weight) reductions on the order of 35%. Table III gives 


information describing the augmented tubes used in this study. 


Table III. Characteristics of Some Commercially-Available 
Augmented Tubes* 


Tube Groove Thermal Enhancement Head 
Geometry Factors Loss 

2 

Factoz 

PiEech Depth Internal External 
Dia Dia 
KG Q 5543 OrsO25 Die lna Lipa OS) 2.03 
KM Oe 54S 0.040 Zook 1.14 3226 
Yalan 05250 (65,0) aba} 1.44 35 1516 
WZ 0.250 CaaS e344 1.94 Le 47 
YM1 Cbs) Os a5 Pet ©) ase 6.02 
YM2 Onee25 OSCO2LS 2.710 1.94 5.69 
Notes: 


1. Tubes KL and KM are products of Wolverine Division, 
UOP Inc. [25]. “Tubes YL, YL2, YMi,“and™ yz) are 
products of Yorkshire Imperial Metals, LTD [26]. 


2. Figures shown are ratios of head loss values to 
the smooth-tube case for the baseline conditions. 
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Figure 8. Effects’of internal and external enhancement on 
tube length; paseline condenser. 
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The smooth-tube baseline is, as is typical, controlled 
in thermal resistance by the tube side film coefficient. For 
this reason internal enhancement is relatively more effective 
(about 2:1) than external enhancement. Figure 8b is a cross- 
plot of Fig. 8a and indicates the combinations of internal 
and external enhancement necessary for a given length reduc- 
tion. Here again it is seen that internal enhancement is 
justifiable to large values relative to external enhancement 
or, by way of example, external enhancement beyond 1.5 may 
not be warranted unless the internal film coefficient is en- 


hanced by a factor of at least 2. 


Variations in Thermal Resistance through the Bundle 


The importance of considering the limiting thermal re- 
sistance is further illustrated in Figs. 9a and 9b. Figure 
9a illustrates the row-by-row variation in internal and ex- 
ternal resistances for a condenser with no external enhance- 
ment. Without internal enhancement, the thermal circuit is 
unbalanced (dominated by the internal resistance) throughout 
the condenser. External enhancement in this case would clearly 
be gilding the thermal lily. With an internal enhancement of 
2.0, however, the two film resistances are relatively balanced. 
Such would be the approximate case for the KL tube of Fig. 8. 

If external enhancement is to be applied, on the other 


hand, such as in the case of YM1l in Fig. 8, additional interna 


46 


Aaqunu MOY 


*soouezSTSerT Tewr3uy 
JO suotjzetTzeA MOor-Aq-Moy. 


°6 oanhty 


Jaquinu Moy 


OL 79 UL 


02 


N4d/do 24 4H 2011 “Yy 


SS GC Cm 


Ste 
BNI 
- 
De 


- —_ => P= 


“Tee — TT 407995 


0°S 


.e) 0°9 
(T= a) 
qUusWwsoUPYyUS 
TeurIsazXS ON (ke) 


Or £ 


0°8 


N3G/4o 534 AU, OL “Y ‘@oURZSLSe4 | eudayl 


47 


dihaneenene is to be recommended. Such a case is illustrated 
in Fig. 9b where, again, the relative balance is obtained 
between the two film resistances by using an inside enhance- 
ment of E. = 3.0 and an outside enhancement of EG = Ls Sica 
should be noted that in Fig. 9b the film resistances have 
been lowered below the thermal resistances due to fouling 
(R.) and tube wall resistance (R,) (resistance due to the 
non-condensible gas film is negligible in this case). There- 
fore, further improvements would be most noticeable from re- 
GuctLon in Re and/or Ry 

Figures 9a and 9b also illustrate the complicated nature 
of the shellside resistance. The dashed lines indicate the 
distribution of external resistance in sectors I, IV, and VI 
of the condenser, with sector I at the top of the bundle (center 
ray at 15°) and proceeding clockwise in 30° increments to 
sector VI. In sector I a procession radially inward from 
row 1 to row 30 is in the direction of gravity and each suc- 
cessive row is inundated by the row above - thus the continuous 
increase in Ry in sector I.. In sector VI, however, each Zow 
radially inward is subject to less inundation so that a gen- 
eral improvement is realized from the point of view of con- 
densate inundation. In every sector, however, the vapor 
velocity is decreasing in a radially inward direction and 
the beneficial effect of vapor velocity is correspondingly 


less in regions toward the center of the bundle. In the top 
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sector (sector I) condensate inundation and vapor shear ef- 
fects change in the same direction with resistance increasing 
radially inward. In the bottom sector (sector VI) the inunda- 
tion and vapor shear change oppositely and there is a general 
balance. It will be noted that near the bundle exit (entrance 
to the air cooler section) in sector VI the decreases in vapor 
shear effects no longer overcome the changes due to reduced 
inundation and the resistance begins to decrease. Sector IV 
is also illustrated in these figures and is seen to represent 
something of a mean behavior among the sectors. The spread 


between sectors in the external resistance, R is seen to be 


Oo’ 
reduced in Fig. 9b because of reduction in the relative effect 


of vapor shear with external enhancement. 


Influence of Vapor Shear and Condensate Inundation 


Figures 10a and 10b illustrate the potentially deceptive 
nature of the influence of vapor shear and condensate inunda- 
tion in condenser bundles. Figure 10a examines the relative 
influence of the two effects upon required tube length and 
their dependency on internal enhancement. In this case the 
favorable effect of vapor shear is roughly balanced by that 
of inundation and in terms of required bundle length they are 
by no means independent of each other. Without enhancement 
the two effects are essentially equivalent so that the HEI 


method, which neglects them both, appears to be adequate (for 
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rather serendipitous reasons). With enhancement, however, 
when row-by-row thermal performance must be considered, the 
neglect of either vapor shear or condensate inundation can 
lead to erroneous results (e.g., the neglect of vapor shear 
would indicate no improvement with an internal enhancement of 
about 1.7). Figure 10b shows that when external enhancement 
is large, these effects are less important because of the 
dominance of thermal resistances other than those due to the 


shellside film. 
Enhancement in Bundles of Various Diameters 


The number of possible comparisons such as those shown 
here is virtually limitless. As a final example, it might 
be useful to illustrate some of the influences to be ex- 
pected with enhancement when the number of tubes in the 
bundle and hence the bundle diameter is allowed to vary. 
Figures lla-llf illustrate this effect for a range of in- 
ternal enhancement factors with no external enhancement (Figs. 
lla-llc) and with Ey = 3.0) (Bigs. Lild-1v£). ’ En each set o£ 
figures the changes in bundle length, volume, and weight are 
shown for variations of + 20% in the number of tubes in the 
bundle with the cooling water velocity held constant at 8.0 
ft/sec. The reader will draw his own conclusions but a few 


interesting points may be highlighted here. 
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Comparing Figures a and d, b. and e, andc and f, it will 
be seen that bundle length, volume, and weight may all be 
reduced through enhancement (E; ee 3.0)" by factors of 
about 38% with the number of tubes held constant at the base= 
line value of 3292. 

An increase in the number of tubes will require an in-= 
crease in the bundle diameter while allowing, at the same 
time, a decrease in the tube length necessary to sustain the 
given heat load. These two effects are in opposition as re=- 
gards bundle volume and weight and in the case of no enhance- 
ment (Figs. lla-llc, Ez. = 1.0) the decrease in length (about 
27%) is sufficient to allow decreases in both volume (203%) and 
weight (13%). At high levels of internal enhancement, how- 
ever, the allowable length reduction due to additional tubes 
is somewhat less and although volume reduction is still ob- 
tained the trend in bundle weight decreases and is eventually 
reversed with E. =\i2e5- and Ey = 1. For example, the length 
reduction accompanying a 20% increase in tubes is about 11%, 
with a volume reduction of 5.4% and a negligible change in 
weight. At even higher values of Es it is seen that a 20% 
increase in tubes leads to an increase in bundle weight - an 
effect that is opposite to that for the situation without en- 
hancement. Reversals in trends, such as these, are often 
impossible to predict if steam-side interaction effects are 


not taken into account. 


If, on the other hand, it is desirable to reduce the 
number of tubes (for space considerations or to reduce cool- 
ing water flow rate, for example) a 20% decrease in the 
number of tubes will give a 10% reduction in bundle diameter. 
Without enhancement, this would require increases in length, 
volume, and weight of 90%, 70%, and 26%, respectively. If 
enhancement (E. = Eo = 3.0) is utilized, however, the per- 
centage increases in length and volume are only 15% and 6% 
while the bundle weight is decreased by about 7%. Thus, en- 
hancement can allow condenser modifications for "packaging" 
purposes that might otherwise be impossible. 

In the examples cited in Fig. ll, the condenser cooling 
water velocity has been held constant at 8.0 ft/s. For 
variable bundle diameters, the cooling water flow rate must 
therefore change in proportion with the number of tubes. 
These variations and constraints are used here for illustra- 
tive purposes and they represent only a few of the cases 
under consideration. It should also be remarked that here 
we have made only minor mention of the changes in cooling 
water head loss and pumping power that accompany the many 
variations shown in Figs. 8-ll. Obviously, the relative 
Merits of adjustments such as these depend upon the premiums 
to be placed upon the costs and benefits. If increases in 
bundle diameter can be tolerated, for instance, there are 
Significant payoffs to be had in reductions in length, volume, 


weight, and cooling water pumping power. 
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SUMMARY 


The authors wish to again remind the reader that these 
calculations are most trustworthy when viewed as bases for 
comparisons. The extremely complex nature of the £ low and 
heat transfer in condenser bundles leaves much to be accom- 
plished through further research. It is hoped that, with 
these precautions, the utility of a code such as MORCON has 
been illustrated. It is expected that future design com- 
parisons such as these will be expedited by the development 
of a numerical optimization code, based upon MORCON, to al- 
leviate the serendipity required to optimize such complex 
phenomena as condenser bundles. 

In closing, it is to be emphasized that several other 
performance factors and design variables, not mentioned here, 
are under study. The various tube-side and shell-side pres- 
sure drop legacies associated with enhancement are calculated 
and these constitute some of the "costs" to be balanced 
against the benefits of enhancement. A paramount factor 
yet to be understood is the extent to which fouling offsets 
the benefits of enhancement. Other design variables, such 
as tube material, tube diameter and wall thickness, condenser 
steam inlet pressure, and cooling water velocity and flow 
rate are under consideration. A change to titanium tubes 
will, in itself, generate a reduction in bundle weight by 


a factor of approximately 3 because of reduced metal density 
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and tube wall thicknesses. With all design variables left 
to vary within reasonable constraints, it is expected that 
improvements in tube bundle weight and volume may be realized 
in excess of those cited here and with acceptable penalties 


in terms of pressure drops and pumping powers. 


CONCLUDING REMARKS 


On the basis of work reported here and elsewhere, it is 
apparent that the continued advancement of the technology 
associated with surface condensers will require improvements 
in the understanding and definition of the flow conditions on 
the two-phase multi-component side of the problem. Whereas 
the prediction of the performance of single tubes - with 
vapor-side conditions well defined - has seemed adequate 
in the past, this situation may be attributed to conservative 
design attitudes and to the fact that thermal resistances 
other than those on the vapor side have dominated the overall 
thermal circuit. Advanced condenser designs that seek optimal 
performance (minimum weight, maximum efficiency, etc.) will 
require a level of analytical sophistication that exceeds 
present-day standards. 

A case in fact is that in which enhanced tubes are con- 
templated for use in condenser bundles. The prediction of 
condenser performance based upon single-tube theory and ex- 


periment cannot be expected to be reliable. Overall condenser 
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performance may exceed or fall short of such predictions 
depending upon the selection of the tubes themselves as 
well as their geometry and location in the condenser bundle. 
Computer-based analyses can be used to advantage in 
the estimation of the tube-to-tube variation of flow condi- 
tions within a condenser bundle. In this report, a condenser 
computer code of moderate complexity has been described and 
its use has been demonstrated in several sample cases. The 
results of these computer experiments have served to illus- 
trate the importance of internal detail in the prediction of 
the performance of condensers utilizing enhanced tubes. In 
addition, these studies have shown that judicious use of en- 
hanced tubes can lead to considerable improvements in con- 
denser performance as well as the realization of heretofore 


unattainable design options. 
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